Multimode uncertainty relations and separability of continuous variable states 
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A multimode uncertainty relation (generalising the Robertson-Schrodinger relation) is derived as 
a necessary constraint on the second moments of n pairs of canonical operators. In turn, necessary 
conditions for the separability of multimode continuous variable states under (m + n)-mode bipar- 
titions are derived from the uncertainty relation. These conditions are proven to be necessary and 
sufficient for (1 + n)-mode Gaussian states and for (m + n)-mode bisymmetric Gaussian states. 
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Quantum mechanical uncertainty principles stem di- 
rectly from the noncommutativity of quantum obscrv- 
ablcs and from the probabilistic interpretation of the 
wavefunction. In a close sense, uncertainty relations 
IH IE Hi are the 'operational' expressions of such fun- 
damental axiomatic features. Besides this outstanding 
fundamental interest, uncertainty relations have recently 
acquired a technical interest of major relevance for quan- 
tum information science 0]. In fact, for any state g of 
a bipartite quantum system, the positivity of the par- 
tially transposed density matrix g (obtained from g by 
transposing the Hilbcrt space of only one of the two sub- 
systems) is a necessary condition for the state to be sepa- 
rable 0,0]. In other words, the violation of the positivity 
of g is a proof of the presence of quantum entanglement 
in the state g, which is thus suitable for various quantum 
informational aims. Now, because uncertainty relations 
for quantum observables derive from the positivity of the 
generic density matrix g, the partial transposition of any 
uncertainty relation (i.e. the relation formally derived 
from the uncertainty relation by considering the partially 
transposed state g instead of g) provides a sufficient con- 
dition for the state g to be entangled. Often, as is the case 
for the continuous variable systems to which this paper is 
devoted, such conditions come in a simple form in terms 
of observable quantities, and are therefore of great ex- 
perimental relevance Q. Furthermore, these conditions 
turn out to be also necessary for entanglement in several 
instances, namely whenever the positivity of the partial 
transpose is also sufficient for separability 0, 0, S IH3 • 
Notably, the relationship between uncertainty relations 
and entanglement is exhibited in discrete variable sys- 
tems as well (see Refs. 0, 0| and, for 'spin-squeezed' 
states, 

In this paper, canonical systems of many modes (like 
discrete bosonic fields in second quantization or motional 
degrees of freedom of material particles in first quantiza- 
tion) are considered. Quantities invariant under linear 
canonical (i.e. symplectic) operations on the field modes 
will be constructed as functions of the second moments of 
the field operators. In terms of such invariant quantities, 
a general uncertainty relation for the second moments of 



any n-mode system will be derived. The partial transpo- 
sition of such a relation will lead to useful entanglement 
conditions for any bipartition of the modes. Crucially, 
these multimode conditions are both easily checkable and 
straightforwardly related to experimental data. 

Uncertainty relations for canonical systems. Let 
us consider a continuous variable (CV) quantum me- 
chanical system described by n pairs of canonically 
conjugated operators {xj,Pj} with continuous spectra. 
Grouping the canonical operators together in the vector 
R = (xi,pi, . . . , x n ,p n ) T allows to compactly express the 
canonical commutation relations (CCR) as [Rj,Rk] = 
2ifljk ^3, where the symplectic form fl is defined as 
CI = ®"c<j with lo = is y (s y standing for the y Pauli ma- 
trix) . Dynamical evolutions of the system must preserve 
the CCR. In particular, transformations acting linearly 
on the vector of operators R (in Hciscnberg picture) must 
preserve the symplectic form Q under congruence. Such 
transformations form the real symplectic group Sp2n,u : 
S e SL(2n,M) : S <E Sp 2n ,n <^ S T nS = n. These 
transformations can be surjectively mapped onto unitary 
operations generated by second order polynomials in the 
canonical operators (metaplectic representation). 

Any state of an n-mode CV system is described by 
a positive, trace-class operator g. Let us define the 
2n x 2n matrix of second moments, or "covariance ma- 
trix" (CM), er (with entries cr;.j) of the state g as cr;.j = 
Tr [{R i ,R j }g}/2 - Tr [Rig\Tr [Rjpl The full uncertainty 
relation for such a system reads |l4j 



o- + m > 



(i) 



This inequality - which derives solely from the CCR and 
from g > - is the only condition a symmetric 2n x 2n 
matrix has to satisfy to qualify as the bona fide CM of 
a physical state. Because of the skew-symmetry of O, 
Eq. Q ensures the definite positivity of cr: a > 0. For 
future convenience, let us define the 2x2 submatrices 
7^ [with entries ("1ij)h,k] of the CM cr as (j l3 )h,k = 
(T(j + ^_ 1 ) ^j + f t _i) (each of them describing one mode or 
the correlations between one pair of modes). 

For n = 1, the uncertainty principle Q reduces to 
Deter > 1 (i.e. to the "Robertson-Schrodinger" uncer- 
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tainty relation |2j). This relation along with the condi- 
tion er > are equivalent to the uncertainty relation 
for single-mode systems. For a two-mode system, one has 

Dct a + 1 > Ai , (2) 

where Ai = J2i j^^^lij- Note that the quantities 
Deter and Ai, entering in the previous Ine qua lities, are 
invariant under symplectic transformations |l5j . The in- 
variant nature of the uncertainty principle implies that, 
indeed, the expression of the uncertainty relation for a 
general n-mode CM er must be possible in terms of sym- 
plectic invariants constructed from the entries of er. 

Construction of the symplectic invariants. - Let us 
recall that, because of the definite positivity of er, one 
can apply a seminal result by Williamson 0] to the 
quadratic form er to infer the following basic result: for 
any CM er there exists a (non-unique) symplectic trans- 
formation S G Sp2n.B. such that S 1 erS = v , where 
v = diag (i>j, Vj) . The quantities {vj} are referred 
to as symplectic eigenvalues, whereas the matrix u is the 
"normal form" of the CM er. The uncertainty principle 
can be equivalently recast in terms of the {vj} as 

Vj>l for j = 1, . . . , n . (3) 

The symplectic eigenvalues are clearly symplectic in- 
variants but, for an n-mode system, their analytical ex- 
pression in terms of the second moments turns out to be 
rather cumbersome (when possible at all). Indeed, the 
symplectic eigenvalues can be computed by diagonalising 
the matrix fier, whose eigenvalues turn out to be {TWj} 
for j = 1, . . . , n. The latter statement is easily proved 
by checking it on the normal form u and by considering 
that Q,u = ttS T erS = S^ 1 i}crS for some S G Sp 2n ,]R- 
Now, a natural choice of symplectic invariants, dictated 
by the previous equation, is given by the principal minors 
of the matrix fier, manifestly invariant under symplectic 
transformations acting by congruence on er. Let MjJa) 
be the principal minor of order k of the matrix a |17| . 
then the symplectic invariants of an n-mode state {A™} 
for j = 1, ... ,n are defined as 

A" = M 2 j(Oer) . (4) 

The principal minors of odd order vanish because of the 
alternate sign in the spectrum of fier, thus leaving us with 
n independent symplectic invariants {A"} (as many as 
the symplectic eigenvalues). The quantities {A™} are 
also known as "quantum universal invariants" |18| . The 
expression of the invariants {A™} in terms of the sym- 
plectic eigenvalues {vj} can be retrieved by considering 
the normal form u and reads 

a;-EII^ (5) 



where the sum runs over all the possible j-subsets iS™ of 
the first n natural integers (i.e. over all the possible com- 
binations of j integers smaller or equal than n) . Clearly, 
one has A™ = Deter while, for two-mode states, the in- 
variant Ai coincides with the quantity Ai appearing in 
the uncertainty relation . 

Symplectic uncertainty relations. - Let us consider an 
n-mode CV system and define the quantity E ra as 

n 

£„=^(-ir +j 'A?, (6) 

where we assume Aft = 1. Now, Eqs. 10 and imply 

n 

E n = 11(^-1). (7) 

i=i 

Thus, Inequality 10 leads to the following statement: 
Symplectic uncertainty relation. Let er be the covariancc 
matrix of an n-mode continuous variable state. The sym- 
plectic invariant E„, determined according to Eqs. 
and 10, fulfills the inequality 

E„ > . (8) 

Inequality (0, reducing to the well known relations 
Det er > 1 and (0 for, respectively, n = 1 and n = 2, pro- 
vides a general way of expressing a necessary uncertainty 
relation constraining the symplectic invariants. Eq. (0 
shows that, actually, Inequality is only necessary and 
not sufficient for the full uncertainty relation to be 
satisfied, as it is not able to detect unphysical CMs for 
which an even number of symplectic eigenvalues violates 
Inequality 0. This impossibility is not due to any fun- 
damental lack of information in the symplectic invariants, 
since their knowledge allows to determine the symplec- 
tic eigenvalues as the n solutions for v of the following 
system in the unknown {A™ -1 , v\ (for j = 1, . . . , n — 1): 
A; 1 = v 2 A™zl + A™" 1 for j = l,..., n (A^ 1 = 0). For 
ri = 2, the additional proviso Det er > 1 is enough to rule 
out the undetectable case so that the uncertainty prin- 
ciple is equivalent, for two-mode states, to the set of 
conditions A| - Af + Ag > , A| > 1 , er > 0. 

Symplectic separability criteria. - The positivity of 
the partially transposed state ("PPT criterion") is a 
necessary condition for the separability of any bipartite 
quantum state. Conversely, the violation of such posi- 
tivity is a sufficient condition for a quantum state to be 
entangled. Moreover, as far as the CV systems here ad- 
dressed are concerned, the PPT criterion turns out to 
be sufficient as well for the separability of (1 + ?i)-mode 
Gaussian states (i.e. states with Gaussian Wigner func- 
tion) and of bisymmetric (m + n)-mode Gaussian states 
(here and in what follows, we refer to a bipartite '(m+n)- 
mode' CV state as to a state separated into a subsystem 
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A of m modes and a subsystem B of n modes). These 
facts come in especially handy for CV systems, as the 
action of partial transposition on covariance matrices is 
easily described. Let g be a (to + n)-modc bipartite CV 
state with 2(m + n)-dimensional CM cr. Then the CM cr 
of the partially transposed state g with respect to, say, 
subsystem A, is obtained by switching the signs of the 
to momenta {pj} belonging to subsystem A: cr = TcrT 
with T = ©™s 2 © J-2n, where lan and s z stand for the 
2?? x 2n identity matrix and for the z Pauli matrix. Now, 
in analogy with Inequality Q derived from the positivity 
of the density matrix g, a (generally) sufficient condition 
for separability derived by the PPT criterion is given by 
cr + ifl > Sj or, in terms of the symplectic eigenvalues 
{Dj} of the partially transposed CM cr (whose normal 
form will be henceforth denoted by £>), as 

v 3 > 1 . (9) 

The previous findings allow to recast such separabil- 
ity criteria for (to + n)-mode states in terms of par- 
tially transposed symplectic invariants {A.,}, defined by 
Aj = M2j(ficr). This simple result will be precious: 
Little lemma. Let cr be the physical CM of a state of a 
(m + n)-mode CV system, with to < n. Let er = TcrT 
be the partial transposition of cr with respect to any of 
the two subsystems. Then, at most to of the symplectic 
eigenvalues {vj} of er can violate Inequality ©. 

Proof. Suppose the transposition is performed in the 
rn-mode subsystem. Let T>(a) be the dimension of the 
subspacc upon which the generic matrix a is negative 
definite. Since T reduces to the identity on a (2n + m)- 
dimensional subspace, Inequality cr + iQ, > reduces to 
the (definitely satisfied) Inequality QJ on such a sub- 
space, thus implying T>(cr + if2) < to. One has then 
T>(D + iil) = T>(cr + ifl) < to, because i> + iil — 
S T (a + i£l)S for S £ Sp 2 (rn+n),R an d the signature is 
preserved under congruence. The eigenvalues of i> + iQ 
are given by {9j =F 1}, thus proving the lemma as the 
{i>j} have to be positive {TcrT > because cr > 0). The 
choice of the transposed subsystem is not relevant, since 
the action by congruence of the matrix ffi™ +n s z turns cr 
into the partial transpose under the n-mode subsystem 
cr' and ft into -f2, and cr' - iVl > cr' + ifl > 0. □ 

In analogy with Eq. ©, let us define the transposed 
invariant t m+n = E^"(-l) m+ " +i A^-. The inequality 

t m+n > , (10) 

being necessary for Inequality to be satisfied, is a 
necessary condition for separability under (to + n)-mode 
bipartitions and is thus a sufficient condition to detect 
entanglement in such a multimode system^ irrespective of 
the nature of the state under examination |l9j . Inequality 
cannot detect the negativity of the partial transpose 
whenever an even number of symplectic eigenvalues vi- 
olates condition 0. However, because of the previous 



lemma, for (1 + n)-modc Gaussian states (for which the 
PPT criterion is necessary and sufficient for separability 
0), at most one partially transposed symplectic eigen- 
value can violate Inequality Q. Inequality (f 1 Of) is then 
always capable of detecting such a violation. The same 
argument applies to 'bisymmetric' Gaussian states, de- 
fined as the (to + n)-mode Gaussian states which are 
invariant under mode permutations internal to the rn- 
mode and n-mode subsystems. A bisymmetric Gaussian 
state with CM cr can be reduced, by local symplectic 
operations (on the to- mode and n-mode subsytems), to 
the tensor product of a two-mode Gaussian state and 
of uncorrelatcd thermal states with global CM <T2: 
(72 = S T crS for some S £ Sp2m,n © Sp2n,n- The lemma 
above implies T>(TcrT + ifl) < 1 from which, observing 
that (TS J T) (TcrT + iO) (TST) = Tcr 2 T + in, one infers 
that at most one partially transposed symplectic eigen- 
value of the CM cr can violate Inequality 10 ■ Notice that 
the locality of the operation S is crucial in establishing 
this result, as it implies (TS T T)fl(TST) = Q. 

Summing up, Inequality $10\) is necessary and suffi- 
cient for the separability of (1 + n)-mode and of bisym- 
metric (to + n)-mode Gaussian states |2l| . It has to be 
noted that an effective strategy - based on an iterative 
ma p - to decide the separability of such states does ex- 
ist |23|. However, the condition l|10|) on the second mo- 
ments can be readily analytically verified and may thus 
be very helpful in the detection of interesting CV entan- 
gled states. This is especially relevant in view of the cur- 
rent experimental developments in the implementation 
of multipartite CV protocols |24|j which heavily rely on 
symmetric resources (perfectly discriminated by the pre- 
vious condition). A concrete example might better em- 
phasize the usefulness of relation (|1C)|) . The "GHZ-type" 
states considered in Refs. 0113, obtained by inserting 
squeezed vacua into an array of beam splitters, are the 
prototypical resource for the implementation of telcpor- 
tation networks. Moreover, they turn out to be the sym- 
metric Gaussian states maximising both the couplewisc 
(between any pair of modes) and the genuine multipar- 
tite entanglement [2(| . To fix ideas, let us consider a four 
mode state [whose CM is given by Eq. (31) of Ref. for 
N = 4] with squeezing parameters r\ = r2 = r (see the 
notation of Ref. 0)- Realistically, thermal noise with 
mean photon number q—1 will also be assumed to affect 
the creation of such states (this amounts to multiplying 
the CM by q). Here, two kinds of bipartitions are of 
interest: the one between two modes and the remain- 
ing two and the one between one mode and the other 
three. Working out the partially transposed symplectic 
eigenvalues to check the separability of such bipartitions 
would require the numerical solution of an eighth degree 
algebraic equation. Instead, Inequality (|10|) allows to an- 
alytically verify the separability of any bipartition with 
the elementary operations (multiplication and addition) 
needed to compute the minors A^. For instance, for 2x2 
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bipartitions one gets Aj = q 2j (1 + g cosh(Ar))(C 4 — 2g), 
where g = min(j, 4— j) and Cj is the binomial coefficient. 
Checking condition (|10|l is then straightforward and leads 
to the following analytical relation for the separability 
of the considered bipartiton: cosh(4r) < (q 4 + l)/(2q 2 ). 
Clearly, the advantage of such a strategy gets more and 
more relevant as the number of modes increases. 

The compact and elegant forms in which uncertainty 
relations and separability criteria (remarkably necessary 
and sufficient for two relevant classes of Gaussian states) 
have been recast in the present paper were obtained re- 
lying exclusively on symplectic analysis applied at the 
pha se sp ace level. Much scope is left to such a longstand- 
ing [l4J kind of approach, in particular concerning the 
entanglement characterization of CV states. The analy- 
sis of Gaussian states based on symplectic invariants has 
provided remarkable insight into the entanglement prop- 
erties of two-mode states [2t| and could be, employing 
the techniques developed in this paper, extended to mul- 
timode settings and to the analysis of multipartite contin- 
uous variable entanglement |7l l24Ll25ll28l. which has been 
lately drawing considerable attention ill Hill HII. At 
a more fundamental level, this paper is intended to shed 
further light on the constraints imposed by quantum me- 
chanics on second moments of canonical operators and 
on the symplectic structure underlying the evolution of 
bosonic modes under quadratic Hamiltonians. 
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